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We calculate for the first time the probability density func- 
tions (PDFs) P of burst energy e, duration T and inter-burst 
interval r for a known turbulent system in nature. Bursts 
in the earth-sun component of the Poynting flux at 1 AU in 
the solar wind were measured using the MFI and SWE ex- 
periments on the NASA WIND spacecraft. We find P(e) and 
P(T) to be power laws, consistent with self-organised critical- 
ity (SOC). We find also a power law form for P(r) that distin- 
guishes this turbulent cascade from the exponential P(r) of 
ideal SOC, but not from some other SOC-like sandpile mod- 
els. We discuss the implications for the relation between SOC 
and turbulence. 



In their seminal papers |],§, Bak et al. (BTW) demon- 
strated that a discrete cellular automaton model of an 
artificial sandpile had a spatial response to slow fuelling 
that was characterised by a scale-free distribution of en- 
ergy release events or "avalanches" (see also ||). Scale 
invariance was shown by a power law probability density 
function (PDF) P of avalanche area A, P{A) = CA~ a . 
This scale-invariant spatial structure led BTW jy to pro- 
pose the sandpile as a toy model of turbulence because, 
in Kolmogorov turbulence jj], long- wavelength, injection- 
range perturbations cause a scale-free forward cascade of 
energy transport until the dissipation scale is reached and 
therefore one might expect the PDFs of burst quantities 
in turbulent systems to be power laws too. These have 
recently been shown in burst area A for a generic inverse 
cascade model H, in burst energy e and duration T for 
both a shell model |6) and reduced 2D MHD turbulence 
simulations J7j , and in peak burst power for ID MHD 
turbulence ||T 

Boffctta et al. f| (hereafter B99) have also shown that 
the PDF P of inter-burst intervals r in a shell model of 
turbulence is a power law too but that this is not so for 
the BTW sandpile in which P(t) is exponential. B99 
postulated that the power law P(r) found for solar flares 
P| was consistent with a shell model of turbulence rather 
than the BTW sandpile. Here we demonstrate for the 
first time that the predicted avalanche phenomenology 
(power laws in P(e),P(T) and P(r)) of a shell model of 
turbulence is observed within a natural system - the solar 
wind - for which there is direct independent evidence of 
turbulence 

The solar wind is a near-radial supersonic plasma out- 
flow from the solar corona which carries with it solar 



magnetic flux into interplanetary space by virtue of the 
plasma's very high electrical conductivity. In this ideal 
magnetohydrodynamic (MHD) approximation, the elec- 
tric field E in the rest frame of the moving plasma is 
given byE = E + vxB = from Ohm's law. The 
electromagnetic energy (Poynting) flux E x H along the 
sun-earth line x can be approximated by v [By + B%) / no 
assuming a radial solar wind. This quantity was calcu- 
lated from "key parameter" measurements of B and | v | 
from the MFI and SWE experiments, respectively, on the 
WIND spacecraft |11 between January 1995 and Decem- 
ber 1998 inclusive. The typically 80-100 s averaged mea- 
surements of | v | were interpolated onto the 46 s time 
samples of B. 

In the resulting time series, bursts were identified, by 
the method used in [jl2|, as intervals when the Poynt- 
ing flux exceeded a given fixed threshold. Thresholds 
were set at the 10,20, ...90 percentiles of the cumulative 
probability distribution of the Poynting flux. For each 
threshold, the PDF of the burst energy e, burst lifetime 
T, and inter-burst interval r was calculated, where the 
burst energy is the sum of the Poynting flux samples over 
the burst lifetime T. The PDFs are shown in Figure 1. 
The burst energy PDF (top panel) can be seen to have 
a power law region over about 4 orders of magnitude be- 
tween about 10 -5 and 10 _1 J m -2 . The burst lifetime 
PDF (second panel) also exhibits a power law region and 
can be fitted by a power law with exponential cut-off sim- 
ilar to that found previously for the solar wind e function 
fUSf . In these respects, the solar wind Poynting flux has 
the avalanche phenomenology common to both the BTW 
sandpile and turbulence. 

The inter-burst interval PDF has been plotted on both 
a log-log scale (third panel) and a log-linear scale (bot- 
tom panel). It is readily seen that this PDF is a power 
law rather than an exponential. A power law with an ex- 
ponent of 1.67 is shown by the thick dashed curve in the 
third panel. This power law form distinguishes the solar 
wind from a system having the properties of the BTW 
sandpile and instead shows it to be consistent with the 
shell model of turbulence used by B99. This is the same 
result they found for solar flares, for which there was not 
the direct independent evidence of turbulence that there 
is for the solar wind. 

It is possible that the solar wind avalanche phe- 
nomenology is simply dominated by the advection of an 
already turbulent fluid from the sun rather than by an 
energy cascade within the solar wind itself (S. C. Chap- 
man, personal communication, 1999). We can expect 
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the solar wind outflow from the sun to be strongly in- 
fluenced by energy dissipation events in the solar corona 
such as nanoflares Q because these events can change 
the thermal pressure gradient that drives the solar wind 
jl5| and/or allow reconfigurations of the solar magnetic 
field that aid or inhibit plasma outflow from the sun. 
These observations are also topical in magnctosphcric 
physics because we have previously shown fl3| a simi- 
larity between the avalanche phenomenology present in 
geomagnetic perturbations (ljjf^] (which measure dissi- 
pative currents in the Earth's ionosphere) and that in 
the energy delivered by the solar wind to the Earth's 
magnetic and plasma environment. Independently, an 
analysis of the R/S Hurst exponents of solar wind vari- 
ables and magnetospheric indices flcfl has drawn similar 
conclusions. 

So what does the observation of avalanche phe- 
nomenology in a natural system tell us about its physics? 
BTW postulated ||,§] that the appearance of "avalanche 
phenomenology" (power law burst PDFs) in Nature 
was due to an underlying fixed point in the dynamics 
( "criticality" ) which was attractive ("self-organised") - 
Self-Organised Criticality (SOC). Renormalisation group 
studies |l7j have demonstrated that the Abelian BTW 
model indeed exhibits an attractive fixed point. However, 
although BTW argued that SOC implied avalanche phe- 
nomenology, the converse is not true; and in particular 
the observation of avalanche phenomenology in natural 
systems JXs) does not by itself prove that such systems 
are SOC. 

There are many examples of systems that are either 
not self-organised or not critical, or both, that neverthe- 
less present avalanche phenomenology. Avalanche phe- 
nomenology has been seen in the forest fire model ]l9| ] 
controlled by a repulsive rather than an attractive fixed 
point; it thus has to be tuned to exhibit scaling p0| . 
Some other models fl8| exhibit power-law distributions 
without finite size scaling and so are not bona fide crit- 
ical. Avalanche phenomenology can also be produced 
by coherent noise driving pl| ] or by "sweeping of an in- 
stability" p^j. In addition, the fixed-threshold method 
of estimating burst sizes that was used in |l2| and the 
present work may generally result in scale-free PDFs if 
applied to certain types of time series. The action of 
slicing through a fractional Brownian motion (fBm) time 
series at a fixed level generates a set of crossing times 
known as an isoset, for which the PDF of the time in- 
terval between two subsequent crossings has a power law 
form p3fl . Hence the burst duration and inter-burst in- 
terval statistics drawn from such an fBm time series by 
the fixed threshold method would also be expected to be 
power laws. 

Clearly it is not sensible to apply the SOC label gener- 
ally to systems exhibiting avalanche phenomenology [p4[ . 
Instead we should follow B99 in using a restricted defi- 
nition of SOC, implicit in BTWs choice of name, as be- 
ing the mechanism of self-organisation to a critical state. 
From this point of view, in order to show the presence 



of SOC, one has to demonstrate those properties of self 
organization and criticality that are unique to the pro- 
cess of SOC rather than simply observing the avalanche 
phenomena that SOC was designed to account for. 

In consequence, the important question remains Q| 
as to the generality of B99's identification of an expo- 
nential P(t) with the SOC mechanism. Exponential 
P(t) implies that energy release episodes are uncorre- 
cted in time because of the standard result that Poisson- 
distributed random numbers have an exponential distri- 
bution of waiting times. This will give rise to a l// 2 
power spectrum [ fL8| for frequencies higher than those 
corresponding to the longest correlation time. In the 
BTW model, this is the time for the longest avalanche 
and is set by the system length. Jensen et al. (2^] found 
that the BTW system had a l// 2 high frequency power 
spectrum in energy flow down the sandpile, rather than 
the 1// spectrum indicative of long-time correlation. 

Whilst exponential P(t) certainly holds for the BTW 
sandpile p6p|], this is not true for some other sand- 
pile models. For example, let us consider the nearest 
neighbour OFC model p7|,|l8|]. The conservative form 



of this model has been shown to be critical |28| and to 
evolve to a steady state 18 . In this case, P(t) is found 
to be exponential |27|. However, there is also a non- 
conservative form of the nearest neighbour OFC model 
p7| , ^8t in which dissipation is introduced. This was re- 
cently shown to cease to be critical |[28| and, in this dis- 
sipative case, P(t) is found to differ from an exponential 
p7| , p9| . This supports the identification of exponential 
P{t) with SOC. 

Three classes of sandpile model, all of which modify 
aspects of BTW SOC, exhibit time correlation between 
bursts - variously reported as a non-exponential P(j) in 
the dissipative OFC model J27]]2|] and as a "1/f power 
spectrum in both running |30| and continuous (e.g. [[h]) 
sandpiles. However, it has yet to be shown that any of 
these systems are still SOC in the sense of both posessing 
an attractive fixed point and showing finite size scaling. 

If B99 are correct in identifying time correlation of 
bursts as a diagnostic for the absence of SOC, then there 
should then be no instance of a model that has an attrac- 
tive fixed point and finite size scaling (self-organized and 
critical) and which also has time correlated bursts of en- 
ergy flow (specifically a 1// spectrum or nonexponential 
P(t)). That is, the time correlation in dissipative, run- 
ning and continuous sandpiles is actually the signature 
of the breakdown of self-organized criticality. The appar- 
ent paradox of the observation of scale-free burst PDFs 
in such models is resolved when one recognises that scal- 
ing may survive away from the fixed point, and can thus 
co-exist with time correlation |2^] . Scaling in both space 
and time can thus be a robust "generic" property of such 
"near-SOC" systems even if exact criticality is not. Bof- 
fetta et al.'s test can then test for the presence of SOC 
but cannot distinguish any of the modified sandpiles from 
turbulence models, and hence such "near-SOC" models 
remain possible descriptions of turbulence. 
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FIG. 1. Probability density functions of burst measures for 
the solar wind Poynting flux. From top to bottom, the mea- 
sures are burst energy e, duration T, and inter-burst interval 
r. The PDFs of all measures have power law regions. 
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